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Abstract. We address the phase structure of color superconducting quark matter at high quark density. 
Under the electric and color neutrality conditions there appear various phases as a result of the Fermi 
surface mismatch among different quark flavors induced by finite strange quark mass; the color-flavor 
locked (CFL) phase where quarks are all energy gapped, the -u-quark superconducting (uSC) phase where 
it-quarks are paired with either d- or s-quarks, the d-quark superconducting (dSC) phase that is the d- 
quark analogue of the uSC phase, the two-flavor superconducting (2SC) phase where it- and d-quarks 
are paired, and the unpaired quark matter (UQM) that is normal quark matter without pairing. Besides 
these possibilities, when the Fermi surface mismatch is large enough to surpass the gap energy, the gapless 
superconducting phases are expected. We focus our discussion on the chromomagnetic instability problem 
related to the gapless CFL (gCFL) onset and explore the instability regions on the phase diagram as a 
function of the temperature and the quark chemical potential. We sketch how to reach stable physical 
states inside the instability regions. 

PACS. 12.38.-t Quantum chromodynamics - 12. 38. Aw General properties of QCD 



1 Family of color superconducting phases 

The phase structure of matter composed of quarks and 
gluons described by Quantum Chromodynamics (QCD) 
has been investigated for many years, and in the high tem- 
perature and low baryon (or quark if deconfined) density 
region which is accessible in heavy-ion collisions interest- 
ing discoveries have been reported both in theories and 
in experiments. In the high density and low temperature 
region, on the other hand, our knowledge is still poor as 
compared with the rich physics expected in this region. 
Heavy-ion collisions are not suitable for the purpose to 
probe dense and cold quark matter, and such a system 
could be realized, if any, only in the cores of compact stel- 
lar objects. The experimental data from the universe is, 
however, quite limited and there is no smoking-gun for 
color superconductivity so far. Nevertheless, the theoret- 
ical challenge to explore the QCD phase structure is of 
great interest on its own. Also it would be potentially im- 
portant in studying the structure and evolution of neutron 
stars. 

In order to look into a dense quark system, some of 
concepts known in condensed matter physics have been 
imported into QCD in hope of analogous phenomena tak- 
ing place. In this sense the physics of dense quark matter 
is, so to speak, "condensed matter physics of QCD" as 
articulated clearly in the review pQ. Superconductivity is 
definitely one of them. In general the Cooper instability 
inevitably occurs wherever there are a sharp Fermi surface 



below which particles are degenerated and an attractive 
interaction between particles on the Fermi surface. Even 
QCD matter is not an exception and the condensation of 
quark Cooper pairs leads to color superconductivity. 

A major difference between ordinary electric supercon- 
ductivity in metals and color superconductivity in quark 
matter arises from the fact that quarks have three colors 
and three flavors in addition to spin, so that quark mat- 
ter allows for many pairing patterns. The color and flavor 
degrees of freedom make the dense QCD phase structure 
so complicated that subtleties still remain veiled. In this 
article we shall argue what has been clarified by now and 
what should be solved in the future mainly following the 
discussions in my recent papers 

1.1 Pairing patterns 

Many theoretical works have revealed that the predomi- 
nant pairing pattern is anti-symmetric in spin (spin zero), 
anti-symmetric in color (color triplet), and anti-symmetric 
in flavor (flavor triplet). Moreover, the color-flavor lock- 
ing is known to be favored in energy, so that the color and 
flavor indices are locked together. Then there are three 
independent diquark condensates or gap parameters fj]; 

(MC-fttf) ~ ^ie abl e m + A 2 e abs e l]2 + A 3 e ab3 e ll3 , (1) 

where {i,j) and (a, b) represent the flavor indices (u, d, s) 
and the color triplet indices (red, green, blue) respectively. 
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Gap Parameters 


Phase 


Ai ± 0, A 2 + 0, A 3 ^ 


CFL 


A x = 0, A 2 0, A 3 yi 


uSC 


Ax j= 0, zA 2 = 0, A a j= 


dSC 


zii = Zi 2 = 0, A :i 


2SC 


/Ai = Ai = A 3 = 


UQM 



Table 1. A family of color superconducting phases; the color- 
flavor locked (CFL) phase, the -u-quark superconducting (uSC) 
phase, the d-quark superconducting (dSC) phase, and the two- 
flavor superconducting (2SC) phase. UQM is unpaired quark 
matter. The sSC, 2SCsu, and 2SCds phases would not appear 
in the QCD phase diagram. 



The charge conjugation C and the Dirac matrix 75 are 
required to make Q a Lorenz scalar. Of course you can 
consider other kind of pairing between quarks which are 
totally anti-symmetric under exchange, and even different 
types of condensates may coexist. Actually diquark con- 
densates such as spin-zero pairing in the color symmet- 
ric (color sextet) channel and spin-one pairing between 
quarks of the same flavor have been analyzed quantita- 
tively [SBHIEIIEI an d known to be much smaller than the 
predominant condensate. In this article we shall simply 
neglect them. 

Under the pairing ansatz l(Tjl. A\ is a gap parameter 
for the Cooper pairing between d and s flavors and green 
and blue colors. That is, A\ is for bd-gs and gd-bs quarks 
and Ai and A-$ are to be understood likewise; 



A 2 
A 3 



bd-gs gd-bs 
rs-bu bs-ru 
gu-rd ru-gd 



Each gap parameter is either zero or finite and there 
are 2 3 = 8 combinations accordingly. Only five of eight 
phase possibilities as listed in Table are of our inter- 
est relevant to the QCD phase diagram. When three gap 
parameters are all nonzero, this state is called the color- 
flavor locked (CFL) phase. When only A\ is zero, this 
is the u-quark superconducting (uSC) phase named after 
the fact that remaining Ai and A% are gap parameters for 
pairing involving u-quarks. The ci-quark superconducting 
(dSC) phase is understood in the same way. The two-flavor 
superconducting (2SC) phase has only A 3 which is non- 
vanishing. The question is; where and how they show up 
on the phase diagram. The next section is devoted to this 
issue. 

It is worth mentioning that these phases can be char- 
acterized by global symmetry breaking patterns. In partic- 
ular the second-order phase transitions between the CFL 
phase and the uSC or dSC phase belong to the universality 
class same as an 0(2) vector model 9 . In QCD, neither 
the sSC (s-quark superconducting), 2SCsu (2SC of s- and 
u-quarks), nor 2SCds (2SC of d- and s-quarks) phase is 
realized actually because any pairing containing massive 
s-quarks is disfavored by the Fermi surface mismatch en- 
ergy. (The solution branch of the gap equations belonging 
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Fig. 1. The energy dispersion relations near the average Fermi 
surface p, for 1 (solid lines) and 2 (dotted lines) species without 
pairing which have a Fermi surface mismatch Sfi. The dashed 
curves represent the energy dispersion relations with pairing 
between 1 and 2 species with the gap parameter A > 8^/2. 
When A < S/J./2 the dispersion relations cross zeros as shown 
by the solid curves and gapless quarks appear inside the block- 
ing region. 



to the 2SCsu phase has been examined in Rcfs. 
and confirmed to cost a larger energy.) 



1.2 Gapless superconductors 

The classification of superconducting phases under the 
pairing ansatz Q is not complete until we will take ac- 
count of the gapless superconducting states. They do not 
break any new symmetry and the phase transition to a 
gapless superconductor exists only at zero temperature |2] • 
In this sense the gapless CFL (gCFL) phase for instance 
which we will closely discuss later is not a totally new 
phase but can be regarded as a variant of the CFL phase 
augmented by the presence of gapless quarks. Before ad- 
dressing the CFL problem, we shall see the simplest ex- 
ample that is actually enough for abstracting the essence. 

Let us assume that there are two species of particles, 
1 and 2, which have the Fermi momenta p, — Sfi/2 and 
/2 + 8fx/2 respectively and they form a Cooper pair. The 
energy dispersion relations without pairing are shown by 
the solid and dotted lines in Fig.Q] In the presence of the 
1-2 pair condensate A, the level repulsion by the energy A 
results in the dispersion relations which are smoothly con- 
nected between the hole state of 2 (or 1) and the particle 
state of 1 (or 2). In this simple model of superconductivity, 
therefore, the quasi-particle energy is expressed as 



e±(p)= y/(p-ji)*+A*±%5n 



(2) 



As is obvious from Fig. ^ as well as from the expres- 
sion J2J) , the dispersion relation comes to cross zeros when 
\^Sfi\ > A. The momentum region from one zero to an- 
other zero of the dispersion relation is called the block- 
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ing momentum region because the pairing within this re- 
gion is hindered by degenerated particles. (In the case of 
Fig- particle-2 is degenerated and particle-1 is absent 
in the blocking region.) Generally once a superconductor 
enters the gapless state, the gap parameter significantly 
decreases with increasing blocking momentum region. 

In the language of physics, the condition for the gap- 
less onset \^Sfj,\ > A means that a Cooper pair is not 
stable energetically. The pairing energy 2A is needed for 
breaking a Cooper pair into two particles and at the same 
time the mismatch energy Sfj, is released by doing so. If the 
mismatch is more expensive than the pairing, two particles 
would no longer form a Cooper pair. Roughly speaking, 
the realization of gapless superconductivity can be seen 
as a weak instability disrupting the Cooper pair only in a 
limited blocking region but not ruining the whole super- 
conductivity. 

We will end this subsection with one more remark. In 
this article (and in some literatures) <5/i is frequently called 
the Fermi surface mismatch. It should be kept in mind 
that this mismatch is for the energy dispersion relations 
without pairing. In the presence of pairing, as observed in 
Fig-ffl the (approximate) Fermi surface is provided by the 
average p,. In other words two energy dispersions e^(p) in 

have a common /2. 



1.3 Strange quark mass effect and neutrality 

We have so far illustrated what the gapless superconduct- 
ing state is like. This gapless phase is often called the 
Sarma phase and was first demonstrated as a metastable 
state by Sarma ^2). The Sarma phase has recently at- 
tracted much interest since it was pointed out that it can 
be stable under the constraint to fix the particle num- 
ber dUIlj. Also in the context of QCD the gapless 2SC 
(g2SC) phase first \T^\ and the gCFL phase next ^B] have 
been shown to be stable under the electric and color neu- 
trality conditions, that is; a potential maximum with re- 
spect to the gap parameters turns to be a minimum if ad- 
ditional gap parameter dependence through the neutrality 
conditions is properly taken into account. The gapless su- 
perconducting state stabilized in this way, however, turns 
out to be unstable after all, which is actually the main 
subject we shall argue spending the last half of this arti- 
cle. 

Anyway, postponing the instability problem for a while, 
let us go on our discussion on the gCFL phase. Now we 
have to understand where dfj, could come from. There are 
two distinct contributions from one origin; one is the di- 
rect M s effect on the energy dispersion relations and the 
other is the induced M s effect through the neutrality con- 
straints. 

The direct effect is simple. In the vicinity of the Fermi 
surface (p ~ /i), the quark energy dispersion relation can 
be well approximated as 

Vp r +Ml-V^p-(v-if-)- (3) 



Thus the strange quark mass effect can be incorporated 
by a chemical potential shift proportional to M 2 /fi. It 
is actually the case that the physics should change as a 
function of M^/fi alone. 

An explanation of the induced M s effect needs some 
knowledge on the neutrality conditions. As far as a chunk 
of quark matter, like the cores of neutron stars which 
are of kilometer size, is concerned, neutrality with re- 
spect to charge associated with any gauge field must be 
required jl7U18U19j . Besides, in the case of color charge in 
QCD, bulk quark matter must be a color singlet as a whole 
system. This is a more stringent condition than neutral- 
ity. The singletness is in principle to be implemented by an 
appropriate projection operator onto neutral states. There 
are many quarks in a macroscopic system and there are 
many thermodynamically equivalent states. If the num- 
ber of singlet states grows in the large volume limit in the 
same way as the number of neutral states, the singletness 
condition can be simply ignored in thermodynamic prop- 
erties. This is in fact the case in color superconducting 
quark matter |20) . Thus only the neutrality conditions for 
electric and color charges suffice for our purpose to explore 
the phase structure. 

In the grand canonical ensemble, the number density of 
particles is specified by a chemical potential. We should 
introduce one chemical potential /x e for the electromag- 
netic charge and eight chemical potentials [L a for the color 
charges which are not gauge- invariant in non-abelian gauge 
theories. For the purpose of imposing neutrality it is enough 
to constrain only two eigenvalues of the color charge since 
zero charge remains to be zero charge if it is rotated by 
any gauge transformations. We thus only have to solve 
the color neutrality conditions with respect to two chem- 
ical potentials /13 and |18j . 

Quark matter with three flavors would be automati- 
cally electric and color neutral if M s is zero, meaning that 
fi e , /i3, and should be of order M 2 //i just like the direct 
M s effect. 

In summary, the direct and induced M s effects are con- 
cisely expressed in a form of the effective chemical poten- 
tials for respective quarks with color a and flavor i as 

Mai = M-Me(Q) ii +/i3(T 3 ) aa 

+ M8 (|T 8 ) aQ -(M 2 ) i4 /2^, (4) 

where Q = diag(|, -|, -±) and M = diag(0, 0, M s ) in 
flavor (u,d,s) space and T3 = diag(|,— ^,0) and -^Tg = 

diag(|,^,— |) in color (red, green, blue) space. From 
the Fermi surface mismatch for all the pairings is easily 
inferred; for the bd-gs pairing for instance the mismatch 
is 

Sf-bd-gs = Pbd — Mgs = 5^3 — M8 + -^t ■ (5) 

In the CFL phase at zero temperature a model-independent 
argument JS] yields fx e = (13 = and = —M s 2 /2fi, and 
together with the above expression, the gapless onset con- 
dition (where gapless quarks begin appearing) is 

Spbd-g. = Mf/fi > 2A X . (6) 



4 



Kenji Fukushima: Phase Structure and Instability Problem in Color Superconductivity 



Numerical calculations in a model study |2*lllf)| have con- 
firmed this onset condition being a good estimate. 

2 Phase diagram of dense quark matter 

We will present the phase diagram first and then look 
closely at each phase in turn. Figure is the phase dia- 
gram obtained in the Nambu-Jona-Lasinio (NJL) model. 
The model parameters are chosen as to yield A = 40 MeV 
at fi = 500 MeV and M s = T = 0. The strange quark 
mass is fixed at M s — 150 MeV which is the lowest esti- 
mate for the M s effect in the intermediate density region. 
In Ref. [21] the authors solved the gap equations for the 
chiral condensates as well as the gap parameters and the 
chemical potentials. In the case of weak coupling, A is 
smaller and thus the critical value of the strange quark 
mass (or chemical potential) on the gCFL onset is smaller 
(or larger), and the gCFL phase is not substantially af- 
fected by the chiral dynamics then. 

Even if the chiral condensate is not dealt with dy- 
namically, a weird structure including a first-order phase 
boundary emerges on the phase diagram for large M^f/J,. 
As discussed in Ref. [2], however, the complicated struc- 
ture seen at large M^/fi is model and parameter depen- 
dent. Unfortunately there is no guiding principle to avoid 
unphysical artifacts in a certain parameter choice and one 
has to try and show all the cases for conclusive analy- 
ses In this article we dare not to touch that subtle 
issue and actually this is the reason why we present the 
phase diagram only for \i > 390 MeV; a first-order phase 
boundary between the gCFL phase and the UQM is found 
at n ~ 388 MeV. 

The phase structure depicted in Fig.^is robust at least 
in a sense that the gross features would not be amended by 
assumptions we have made to draw this figure. Besides the 
chiral dynamics the important physics we have dropped 
is the K °-condensation that could delay the gapless on- 
set [22] and the 't Hooft interaction term that embodies 
the instanton induced Ua(1) breaking interaction in the 
NJL model Lagrangian. These two missing effects are ac- 
tually closely linked. They could change the quantitative 
details of the phase diagram but no qualitative essence 
would be altered. We would not refer to the mixed phase 
possibility here, for it has already been settled [TU] . 

We can see that in the lower density region near the 
gCFL phase the uSC region opens at finite temperature 
while the dSC phase is relevant at higher density. We 
can understand this behavior in a model-independent way 
and the existence of the "doubly critical point" is con- 
cluded [2]- In the subsequent subsections we shall discuss 
the important nature of the gCFL phase and how the 
uSC and dSC phases come out from the Fermi surface 
mismatch, and finally comment upon the doubly critical 
point. 

2.1 gCFL phase 

The gapless onset is located at u ~ 455 MeV where A\ — 
24.9 MeV and /is = —24.7 MeV and one can readily con- 



firm that ti 8 ~ -M^/2u and Mf / ' \i ~ 2A\ as indicated by 
the condition © . The gCFL phase has a definite meaning 
only at zero temperature because there is no clear distinc- 
tion between gapless quarks and thermally excited quarks 
at finite temperature. 

The onset condition © is only for bd-gs quarks with 
Ai but gapless quarks can be present in another quark 
sector as well. In fact the energy dispersion relations of 
rs-bu quarks turn out to be identical to those of bd-gs 
quarks as long as the system remains to be in the CFL 
(not gCFL) phase. Thus rs-bu quarks become gapless at 
M 2 J\i ~ 2A 2 where A 2 = A x holds in the CFL phase. For 
larger M s or smaller a, the blocking momentum region in 
the bd-gs sector becomes wider, while the blocking region 
in the rs-bu sector is severely constrained by the neutrality 
conditions so that it must remain tiny (see Refs. pil3lll6) 
for details). As a result the rs-bu dispersion relations are 
kept to be almost quadratic in the entire gCFL region. 
Therefore only A\ in the bd-gs quark sector significantly 
drops in the gCFL region due to the spreading blocking 
region with increasing M 2 //i. 

Because this is an essential point in considering the 
instability problem later, we shall reiterate here; 

bd-gs pairing with A\ — gapless quarks appear at M 2 //i = 
1A\\ the blocking region increases for larger M^/fj,. 

rs-bu pairing with A% — gapless quarks appear at M 2 / fi = 
2A 2 ; the dispersion relations are kept to be almost quadratic 
in the entire gapless phase. 

gu-rd pairing with A3 — gapped quarks only. 

ru-gd-bs pairing with all Z\'s — gapped quarks only. 

2.2 uSC phase 

The uSC phase results from the presence of the gCFL 
phase at zero temperature. In the gCFL phase, as ex- 
plained in the last section, A\ is significantly reduced than 
A 2 and A3. The pairing between u-d quarks is not re- 
duced but A3 is enhanced in the gCFL side as A\ and 
A 2 decreases [21 ■ (Note that A3 evaluated in the 2SC 
state is in general larger than A3 in the CFL state for 
the same diquark interaction.) In this way the ordering 
A\ < A 2 < A3 is realized. 

It is a well-known fact that the critical temperature 
is of order of the gap parameter at zero temperature. We 
can anticipate that A\ would melt first at finite temper- 
ature right above the gCFL region on the phase diagram. 
This means that the uSC phase where only A\ is zero is 
expected when gCFL matter is heated. Figure [31 clearly 
shows that the phase structure is certainly as anticipated. 

2.3 dSC phase 

It is not the Fermi surface mismatch but the average Fermi 
momentum that is more relevant away from the gCFL re- 
gion (see the expression (|2J) or Fig.flJ. At zero temperature 
the average Fermi momenta are common in all the quark 
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Quark Chemical Potential [MeV] 

Fig. 2. Phase diagram of dense quark matter obtained in the NJL model with the parameters chosen to give A — 40 MeV at 
fi — 500 MeV and T = M s = 0. The strange quark mass M s is fixed at 150 MeV under the assumption that the first-order 
chiral phase transition occurred at some chemical potential below 390 MeV and (ss) oc M a above that would not induce large 
corrections to M s . The doubly critical point and the gCFL onset are marked by a cross and a triangle respectively. 



sectors, that leads to equality in the number of nine (three 
colors and three flavors) quarks, and thus the electric and 
color neutrality is enforced |2.'i|. 



The enforced neutrality at zero temperature is broken 
at small temperatures of a few MeV |2|. At higher tem- 
perature, especially in the vicinity of the critical temper- 
atures, it is a good approximation to estimate /z's in the 
normal phase, i.e., (i e = — Mf/4/i and /13 — fig — 0. Then 
the ordering of the Fermi momenta, /i s < fi u < fid is con- 
cluded, which can be understood in an intuitive way; the 
number of s-quarks is suppressed by M s and so <i-quarks 
should be more abundant than s-quarks to maintain elec- 
tric neutrality. From this, the average Fermi momenta 
should obey the following ordering; p, su < fids < find- 



The gap equation to determine Z\'s contains the mo- 
mentum integration around the Fermi surface which effec- 
tively picks up the density of states at the Fermi momen- 
tum. The larger the density of states is, the greater the gap 
parameter becomes. In this way the ordering A<i < A\ < 
A3 is realized from the average Fermi momenta ordering, 
which follows the presence of the dSC phase accordingly, 
as first discussed in Ref. |24j . 



2.4 Doubly critical point 

The phase boundary on which A\ goes to zero crosses the 
phase boundary on which Ai goes to zero at the "doubly 
critical point" where two phase transitions with respect 
to A\ and A% take place simultaneously. Since the exis- 
tence of the uSC and dSC phases are robust and model- 
independent, so is the doubly critical point. 

We would shortly comment upon a puzzling question 
concerning the doubly critical point. The question is the 
following; what are the effects of gauge field fluctuations 
on the doubly critical point? 

In weak coupling the gauge field fluctuations bring 
about an induced first-order phase transition and the criti- 
cal temperature is shifted j25U26| . Between the CFL phase 
and the uSC or dSC phase, no manifest effects would 
be expected because eight gluons are all massive in both 
phases. Therefore we can conclude that the phase transi- 
tions from the CFL phase toward either the uSC or dSC 
phase is surely of second order belonging to the same uni- 
versality class as an 0(2) vector model. 

The gauge field fluctuations play an important role, on 
the other hand, between the 2SC phase and the uSC or 
dSC phase and the phase transition is forced to be of first 
order. 

The doubly critical point is the point at which two 
phase transitions meet and three gluons become massless 
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around it since it faces the 2SC phase. This suggests that 
the phase transition must be of first order and the critical 
temperature should be shifted at the doubly critical point. 
But how is it possible if the phase boundaries facing the 
CFL phase are not shifted at all? Or, are they shifted near 
the doubly critical point by nearly massless gluons? To 
answer this question, further clarification on the treatment 
of the gauge field fluctuations is waited. 



can contribute to II and U22 ■ Solid lines represent the 
quark propagator for rs and gs quarks in this specific case. 
From this kind of diagrammatic deliberation, it is easy to 
confirm that there is no mixing between gluons except 
A3-Ag mixing. Moreover, some algebra can readily tell us 
that the screening masses are degenerate for the gluon 
pairs (Ai,Az), (^4,^5), and (Aq,Ai). In this article, like 
in Ref. [3], we shall denote A1.2, ^4,5, and Aqj to mean 
either of the degenerated pairs. 



3 Chromomagnetic instability 

The gapless superconductors could be stabilized by the 
neutrality conditions, at least, in the variational space 
spanned by gap parameters and chemical potentials, but 
after all, they have turned out to be unstable in a wider 
space including the gauge fields. The chromomagnetic in- 
stability signifies that the Meissner screening mass squared 
becomes negative, that is, the Meissner mass is imaginary. 
Before addressing the chromomagnetic instability problem 
first pointed out by Huang and Shovkovy [27J in the case 
of the g2SC phase and analyzed in Refs. 3,28,23j for the 
gCFL phase, we shall make a brief overview on the Debye 
and Meissner screening masses in the CFL phase, which 
would be useful to make a deeper insight into the instabil- 
ity problem. Our results are summarized in Fig.EHand the 
goal of this section is to explain what is going on inside 
the instability regions presented by shadowed regions in 
this figure. 



3.1 Debye and Meissner screening masses 

The Debye and Meissner screening masses are the screen- 
ing masses for the longitudinal and transverse gauge fields 
respectively, which are defined by 



m Z5,a/3 = - lin l n a°f3 = 0, 



m 



M, a p = | lim ,(<% - q t qj)n^JuJ = 0, q), 



(7) 

(8) 



where % — qi/\q\ and II £p is the polarization tensor for 
the gauge fields A£ with the color and Lorenz indices de- 
noted by a and /x. At high enough density the polariza- 
tion tensor is dominated by the quark-loop contributions 
alone and so we shall neglect the gluon and ghost loops 
that would not depend on /i at the one-loop level. 

Let us take one example to elaborate the quark-loop 
contributions for each gluon. The A\ and A<i gluons with 
the color indices labeled according to the Gell-Mann ma- 
trices in color space couple red quarks to green quarks and 
vice versa, while the flavor is not changed at any gluon 
vertices. For example, a diagram graphically shown as 



(Example 1) 




3.2 Unstable channels and mixing with photon 

If there is something funny in connection with gapless 
quarks in the quark- loop polarization like Example 1, it 
is expected to happen when both two quark propagators 
involve gapless energy dispersion relations. At the gapless 
onset, in fact, the quark-loop polarization consisting only 
of bd, gs, rs, and bu quarks is divergent, leading to the di- 
vergent contributions to the Debye and Meissner masses 
for Ai t 2, A3, and A 8 gluons. You can easily check that the 
polarization diagrams relevant to the A^s and Aqj gluon 
channels must have at least one of gu, rd, ru, bd, and bs 
quarks which are gapped. 

It should be noted here about mixing between the A 7 
photon and the A3 and Ag gluons. In the symmetric CFL 
phase with M s — 0, there is no mixing with respect to 
A3, but in the present case with M s ^ and thus fi e 7^ 
in the gCFL phase, there is A3-mixing as well due to the 
isospin symmetry breaking. The mass squared matrix for 
A 7 , A3, and As is a 3 x 3 matrix with nonvanishing off- 
diagonal components generating mixing. We will denote 
the eigenmodes of the mass squared matrix by A 7 , A3, 
and Ag. The rotated photon represented by the A 7 field 
is massless all the way because the CFL and gCFL phases 
preserve a rotated electromagnetic U(l) symmetry. Once 
mixing occurs among A 7 , A3 , and As , it is simply a matter 

of convention which eigenmode should be identified as A3 
or Ag. By this reason, though Figure has the instability 
region with the label A3 and Ag, it does not mean that 
both of two eigenmodes suffer from the instability, but the 
fact is actually that only one of them does. 



3.3 Divergences at the gCFL onset 

The divergences in the Debye and Meissner screening masses 
at the gapless onset derive from the density of states which 
is divergent when the energy dispersion relations take a 
quadratic form; e(p) ~ (p — /2) 2 /2A That is, the density 
of states n(p) is given by 



n(p) = Anp 2 



de(p) 



dp 



(9) 



and obviously n(p) — * 00 when the slope of e(p) is zero at 
P = P~ 

The next question is whether the divergence is pos- 
itive or negative in the Debye and Meissner masses. A 
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naive intuition would be that both are positive, for in the 
CFL phase it is well established that m 2 M — \m 2 D should 
hold and one might well consider that they are 

correlated in a similar way even in the presence of M s ^ 0. 
This expectation is partly true but in an unexpected man- 
ner as we will shortly see below. 

To find the relations between the Debye and Meiss- 
ner masses, it is convenient to split the quark propaga- 
tor into four distinct parts. If the quark mass effect is 
approximated as © near the Fermi surface, the energy 
projection operator divides the propagator into the par- 
ticle part that is a function of p — p, and the antipar- 
ticle part that is a function of p + p.. Moreover, in the 
Nambu-Gor'kov formalism, the quark propagator is a 2 x 
2 matrix and its diagonal component is a normal prop- 
agation of particles, and its off-diagonal component is 
an abnormal propagation mediated by diquark conden- 
sates. There are thus four distinct combinations; diagonal- 
particle, diagonal antiparticlc, off-diagonal-particle, and 
off-diagonal-antiparticle. 

The diagram shown in Example I is one example con- 
sisting of only the diagonal propagators. We can construct 
another example in the same gluon channel composed of 
only the off-diagonal propagators as follows; 




gu -~ - - rd (Example 2) 



where the upper propagator is the off-diagonal component 
proportional to A y s connecting ru and gd quarks, and the 
lower one is the off-diagonal component proportional to 
A3 connecting gu and rd quarks. 

It is important to note that the possible singular be- 
havior originating from gapless quarks can reside only in 
the particle parts because antiparticles would never be 
gapless. The Debye mass squared receives the contribu- 
tions only from particles which can be separated into the 
diagonal part, [m^Jdiag, and the off-diagonal part, [m|,] ff. 
The Meissner mass, on the other hand, has all of the 
particle-particle, particle- antiparticle, and antiparticlc- antiparticle 
contributions. In Ref. 3 interesting relations have been 
found; 

[ m M]diag(pp) = "IKldiag ' ( 10 ) 



where ['Ti^ f ] c jiag(pp) is a part of the Meissner mass squared 
coming from the diagrams with two diagonal particle prop- 
agators (like Example 1) and [w^f] ff(pp) with two off- 
diagonal particle propagators (like Example 2). The rela- 
tion (|10f) might be strange at a glance, for we already know 
that m 2 M = \rn 2 D in the CFL phase at M s = 0. Then, how 
can we retrieve the relation m 2 M = ^rajyl Actually, as long 
as the system stays in the CFL side, the diagonal particle- 
antiparticle contribution [™^/]diag(pa) is twice larger than 
[ m A/]diag(pp) and changes the overall sign in the relation 
between the Debye and Meissner masses. The other con- 
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tributions, [mf^diag^a) , [m 2 M }oS{ pa ), and [m| / ] off(aa) are 
all negligibly small. 

It might be surprising that these relations i|10|) and 
fTT|l are satisfied as they are even in the gCFL phase! The 
diagonal and off-diagonal parts both can generate diver- 
gent contributions at the gCFL onset. As far as the QCD 
problem is concerned, the diagonal contributions are al- 
ways larger than the off-diagonal ones, leading to the De- 
bye mass squared diverging positively and Meissner mass 
squared diverging negatively at the gCFL onset. The neg- 
ative sign in (|10|) cannot be compensated in the gCFL 
phase since, unlike in the CFL phase, [m!,f]diag(pa) is no 
longer a match for [m^ / ] diag ( P p ) . 

To put it in another way, the key relation (|10|) can be 
stated as follows: In the diagonal part, we can say that 
the particle-particle loops tend to induce paramagnetism, 
while diamagnetism stems from the particle-antiparticle 
loops. Usually in the superconducting phase, the diamag- 
netic tendency is greater enough to exhibit the Meiss- 
ner effect. In gapless superconductors, however, antiparti- 
cles are never gapless and only the particle-particle loops 
are abnormally enhanced due to the presence of gapless 
quarks. As a result of that, the diamagnetism gives way 
to the chromomagnetic instability. 



3.4 Away from the gCFL onset 

Figure 01 shows the instability regions for respective glu- 
ons. The difference between the Ai$ behavior and the 
A^-Ag, behavior can be understood from the difference be- 
tween the bd-gs and rs-bu quark dispersion relations. It is 
only the diagram shown in Example 1 that causes singu- 
lar behavior in Ai ^ around the gCFL onset. This diagram 
has two quark propagators; one of gapless gs-quarks and 
the other of quadratic rs-quarks. Since rs quarks are kept 
to be almost quadratic in the entire gCFL region, as we 
have explained before, the Ai^-instability extends over 
the whole gCFL region at small temperatures. The insta- 
bility would not persist into regions at higher tempera- 
ture because the quadratic dispersion relations are easily 
affected by thermal excitations. 

In contrast, the same species of quarks constitute the 
A^-Ag instability region. The instability caused by two 
quadratic quark propagations lies in the entire gCFL re- 
gion like the Ai ^-instability but only at tiny temperatures 
of order eV. Instability boundaries at such low tempera- 
tures are not visible actually on the phase diagram. The 
instability caused by two gapless quark propagations is lo- 
calized near the gCFL onset and spreads toward high tem- 
peratures than the ^^-instability because it has nothing 
to do with quadratic quarks. 

Although we would not go further into details, the 
Meissner screening masses evaluated in the g2SC phase [27] 
indicate that the Ag-instability occurs starting at the g2SC 
onset and there also arises the instability for Aj.5,6.7 ghi- 
ons whose boundary is found not at the g2SC onset but 
inside the 2SC region. In our results the A^^ and A§^ 
instability regions are located at large Fermi surface mis- 



match and enter the uSC phase at higher temperature 
and then the 2SC phase farther. We conjecture that these 
instability regions for A45 and A§j would be linked to 
the instability found in the 2SC phase. The nature of the 
instability with respect to A^ and A§.y is presumably 
different from the instability near the gCFL onset that 
we have seen in great details. In the aim of disclosing the 
QCD phase diagram in the intermediate density region, in 
particular, the A^^ and A§^ instability deserves further 
investigation. 

4 Speculations 

This final section is devoted to sketching some specula- 
tions on how to reach the stable states inside the instabil- 
ity regions on the phase diagram, which has been barely 
succeeded so far. 

There are already some attempts to interpret and re- 
solve the chromomagnetic instability problem (32 P33IIM1 
I35j . The most important among them is, in my opin- 
ion, the observation pointed out by Giannakis and Ren 
in Ref. that the instability with respect to gluons can 
be interpreted as the instability toward a plane- wave crys- 
talline superconducting phase. 

The most straightforward interpretation of the chro- 
momagnetic instability is, of course, spontaneous gener- 
ation of the expectation value for the transverse gluons. 
The gauge field itself is, however, not a physical quantity 
depending on the gauge choice. Actually, what Giannakis 
and Ren realized in Ref. is that the spontaneously 
generated gauge fields can be absorbed in the phase of 
the gap parameters by means of the gauge transformation. 
It should be noted that whether the gauge invariance is 
maintained or not does not matter. The gauge transfor- 
mation in this manipulation simply means the change of 
variables. 

The crystalline phase had been already studied 
before the chromomagnetic instability was discovered. It 
is called the "crystalline" phase because the gap parameter 
takes a form of 

A(x) = \A\e tq - x , (12) 

which breaks the translational and rotational invariance. 
The essential point is that there is another quark basis 
where the phase factor of (|12|l vanishes at the price of the 
vector potential arising in the effective action. Therefore, 
taking care of a vector potential turns out to be equivalent 
with dealing with the gap parameter with a phase factor 
corresponding to the given vector potential. 

Let us rephrase the above mentioned idea in a slightly 
different way. Supposing we performed the derivative ex- 
pansion of the thermodynamic potential ^^[/^(x)] which 
is now gauged with gluons, then we have in general, 

Qa[A{x)} ~ Q [A] 

-K ab {A}Tr[(d t + iA*)A*(x)] a [(d t -iA l )A(x)]\(13) 

where a and b represent the color triplet indices and A 1 = 
A l a T a with the color adjoint index a. The average over x is 
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symbolically implied in Tr. Then, an alternative definition 
of the Meissner screening mass is immediately available 
from this thermodynamic potential as 



1^ 



d 2 Q A 



M, a /3 3 dAidA^ 



A=0 



2 K ab (T a ) ca (T^ bd A c A d . 

(14) 



where A(x) is assumed to be spatial constant. It might 
be instructive to see how this expression works actually. 
In the 2SC phase for example, A a oc S a3 and (T 1 )^ = 
{T 2 )a3 = (T 3 )a3 = 0, and hence one can instantly con- 
clude that the A\, A2, and ^3 gluons are not Meissner 
screened. In the g2SC case the A s gluon is unstable, which 
can be stated by the condition k 33 < because m 2 M 88 ~ 
n 33 {A 3 ) 2 . 

Now we shall assume the crystalline superconducting 
phase with l|12|) . The curvature of the thermodynamic po- 
tential with respect to q is 



1 3 

l—l 



d 2 Q 



K ab A a A b 



(15) 



In the 2SC phase the instability condition for q to grow is 
given by k 33 < again. This argument works well only in 
the (g)2SC case in which the Ag-instability can be identi- 
fied as the instability toward a crystalline superconducting 
phase equivalently. Essentially the above is what has been 
articulated in Refs. [321133] . 

The generalization to the (g)CFL problem is possible 
with a simple extension of the crystalline ansatz Ijl2jl . As 
suggested by Giannakis and Ren, we can consider a colored 
crystalline superconducting phase with the gap parameter 
taking a form of 



A(x) = \A\e lTaq ° 



(16) 



then it is almost obvious that the curvature with respect 
to q a is identical with the Meissner mass squared. In this 
sense, we can identify the instability regions in Fig. |3j with 
the colored crystalline phase with nonvanishing q a . 

What we should do next is now apparent; q a are the 
new variational parameters to be determined so as to min- 
imize the thermodynamic potential. This is a quite tough 
task, however. The number of the new variables is five cor- 
responding to Ai : 2, At,5, A e ^, A3, and A s . To be worse, 
the rotational symmetry is broken by the direction of q a , 
and the momentum angle integration cannot simplify in 
evaluating the thermodynamic potential. So, the numeri- 
cal calculations are too time-consuming to be done. More- 
over it is difficult to achieve an enough accuracy to get 
reliable outputs in the multi-dimensional numerical inte- 
gration. Maybe we have to abandon this apparent strat- 
egy and instead need to invent a wiser simplification that 
would not lose the important physics. 

One possibility would be to go to the higher orders in 
the derivative expansion l|13l> . The calculation would be 
feasible because the expansion coefficients can be evalu- 
ated at vanishing q a and the rotational symmetry is not 
broken then. Although we certainly have a chance to find 



the energy minimum with some q a , there is no guaran- 
tee that the next higher order terms can be adequate to 
stabilize the potential. In any case this kind of calculation 
has yet to be performed. 

Finally, if I am excused to speak of my own opinion, 
all these theoretical efforts will be tested in the lattice 
QCD simulation someday when the sign problem at fi- 
nite density will be solved and the lattice spacing can be 
small enough to describe high density matter. The phase 
diagram like Fig.|3Jis to be confirmed then. 

This article is based on the talk given for the new tal- 
ent sessions at International School of Subnuclear Physics 
43rd Course held at Erice in Italy from Aug. 29 to Sep. 7 
in 2005. I thank all the organizers for stimulating lectures 
and discussions at school. 
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